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PREFACE 

In  the  early  days  of  steam  generation  for  power,  saturated  steam 
was  used  a  great  deal,  but  in  order  to  achieve  higher  thermal  efficiencies 
and  reduce  errosion  of  steam  turbine  blading  the  steam  power  industry 
rapidly  extended  steam  conditions  into  the  superheat  region.  Except  for 
a  small  amount  of  data  recorded  for  boiling  in  vertical  tubes  with  free 
convection  no  extensive  study  had  been  made  of  the  heat  transfer  and 
pressure  drop  characteristics  of  wet  steam. 

VJith  the  advent  of  the  nuclear  reactor  (a  highly  concentrated  heat 
source)  and  the  forced  circulation  boiler,  wet  steam  has  again  become  an 
important  energy  transfer  medium.  Basic  data  for  wet  steam  under  forced 
circulation  has  become  very  important.  This  topic  has  been  investigated 
in  limited  ranges  by  several  students  at  the  U.  S.  Naval  Postgraduate 
School:  in  19^3  by  Fisher  and  King  (3)j  in  1$£U  by  Davis  and  Duacsek  (l)j 
and  in  195>5>  by  Nelson  (12).  Dengler  (2)  measured  heat  transfer  coefficients 
and  pressure  drop  of  wet  steam  at  low  pressure  in  a  vertical  tube  at  the 
llassachusettes  Institute  of  Technology  in  1953. 

The  primary  purpose  of  this  investigation  was  to  cover  the  range 
from  0  to  100/j  moisture,  with  particular  attention  to  the  annular  flow 
regionj  to  measure  heat  transfer  and  pressure  drop  characteristics  and 
to  compare  the  pressure  drop  data  with  the  methods  of  prediction  developed 
by  Ilartinelli  et  al.  at  the  University  of  California  (7),  (9),  and  (10). 

The  experimental  work  of  this  thesis  was  carried  out  at  the  United 
States  Naval  Postgraduate  School  from  January  to  March  1956. 

The  writer  wishes  to  express  his  appreciation  for  the  guidance  and 
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assistance  of  Professor  E.  E.  Drucker  in  conducting  experimental  vrork 
and  writing  this  paper,  and  to  Chief  F.  H.  HLechen,  USN  for  his  aid 
in  making  changes  to  the  experimental  set  up  and  for  boiler  and  condenser 
operation  during  the  tests. 
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CHAPTER  I 
INTRODUCTION 

The  problems  of  two  phase  flow  have  long  been  apparent  in  special 
industrial  applications,  but  only  specific  problems  have  been  investi- 
gated leaving  a  dearth  of  general  data.  The   characteristics  of  t-.;o 
phase  flow  that  are  of  interest  in  this  investigation  are  the  local 
film  he it  transfer  coefficient  and  the  pressure  drop  per  unit  length. 

Specific  problems  of  two  phase  flow  have  been  carried  out  in  the 
refrigeration  and  in  the  chemical  process  industry.  The  observed  values 
of  film  heat  transfer  coefficient  are  much  greater  with  two  phase  flow 
than  with  the  liquid  alone  or  the  vapor  alone.  Verschoor  and  Stemerding 
(13)  observed  film  heat  transfer  coefficients  for  two  phase  flow  up 
to  seven  times  that  observed  with  the  liquid  alone,  in  an  investigation 
with  air  and  water.  Dengler  (2)  observed  values  up  to  fifteen  times 
greater  with  wet  steam  at  low  pressure. 

1  r  e  primary  purpose  of  this  investigation  was  to  deter;  line  the 
values  of  film  heat  transfer  coefficient  over  a  wife  range  of  moisture 
contents.   In  a  similar  investigation  of  Pre on  at  low  temperature  Yoder 
and  Dodge  (lit)  observed  a  maximum  of  the  film  heat  transfer  coefficient 
at  about  Ii0,j  dry  Preon  by  weight. 

A  secondary  purpose  of  this  investigation  was  to  measure  the 
pressure  drop  over  a  wide  range  of  moisture  contents  and  compare  the 
results  with  the  prediction  methods  of  riartinelli,  et  al.  (7)>  (9)*  and 
(11).  Generally  the  experiments  of  I-iartinelli  were  carried  out  iso- 
thermally  in  horizontal  pipes. 


The  primary  variables  in  this  investigation  are:  the  mass  rate 
of  flow  of  the  liquid,  the  velocity  of  the  liquid,  the  viscosity  of 
the  liquid,  the  density  of  the  liquid,  the  mass  rate  of  flow  of  vapor, 
the  velocity  of  the  vapor,  the  viscosity  of  the  vapor,  the  density  of 
the  vapor,  and  the  rate  of  heat  transfered  per  unit  area.  In  two  phase 
flow  with  the  equipment  available  the  velocity  of  the  liquid  or  the 
vapor  could  not  be  determined.  With  annular  flow  the  vapor  flows  through 
the  center  of  the  pipe  at  high  velocity  while  the  liquid  flows  relatively' 
slowly  along  the  sides.  Measurements  of  the  volumetric  fraction  of  liquid 
for  two  phase  flow  in  horizontal  pipes  have  been  made  by  Hartinelli  (7),  (9), 
and  (10)  by  operating  two  quick  closing  valves  to  isolate  the  section, 
followed  by  a  rapid  blow  down  to  remove  the  vapor.  The  liquid  was  then 
drained  and  measured.  To  measure  the  water  clinging  to  the  walls  of 
the  test  section  a  volatile  fluid  was  used  to  rinse  the  section.  This 
was  collected j  the  volatile  liquid  was  distilled  off,  and  the  remaining 
water  was  measured.  Another  complicated  method  was  used  by  Dengler  (2) 
to  measure  the  volumetric  fraction  in  a  vertical  test  section.  Radio- 
active tracer  was  added  in  minute  quantity  to  the  liquid.  A  Geiger- 
1-iueller  counter  was  moved  up  and  down  outside  of  the  test-section,  and 
the  resultant  counting-rate  measurement,  together  with  the  local  weight 
fraction  vapor  and  the  liquid  density,  could  be  related  to  obtain  the 
volumetric  fraction.  The  above  work  of  Hartinelli  and  Dengler  make  it 
possible  to  estimate  the  volumetric  fraction  for  this  experiment. 

The  magnitude  of  the  flow  rates  in  pounds  ner  hour  of  the  liquid 
and  vapor  were  measured  by  flow  meters  after  separation  by  a  centrifix. 
The  viscosity  and  density  were  determined  from  the  steam  tables  as  a 


function  of  saturated  pressure.   To  ''etermine  the  heat  transfer  coef- 
ficient:  thermocouples  were  used  to  determine  wall  temper? ture;  and 
the  saturation  temperature  corresponding  to  tho  local  Pressure  was 
used  for  fluid  temperature.   In  this  irvestigation  pressure  was  hold 
constant  exceot  for  three  runs.   Therefore,  viscosity  arid  density  of 
the  liquid  and  the  vapor  were  constant.  Isothermal  runs  were  m^de  to 
determine  pressure  drop  at  hOO,  700,  and  10. 0  pounds  per  hour.  Heated 
runs  were  made  with  one  heat  rate  (78,300  3TU/Hr/Pt2)  at  U00,  700, 
1000,  and  12^0  pounds  per  hour.   The  limitations  of  the  equipment  were 
investigated  and  these  are  included  in  Appendix  II. 

The  vertical  test  section  included  four  separately  heated  six  inch 
sections  in  series  with  an  inside  diameter  of  \   inch.  The  pressure 
drop  was  determined  from  the  vtlue  measured  across  the  two  center 
sections.   The  heat  transfer  coefficient  was  determined  as  the  average 
heat  transfer  coefficient  for  the  two  center  test  sections. 

The  results  of  this  experiment  were  correlated  with  percent  moisture 

as  well  as  with  X. .  (a  oimensionless  parameter  developed  by  Martinelli 

(7)  ).  Heat  transfer  coefficients  and  pressure  drops  are  usually  related 

in  single  phase  investigations  to  several  oimensionless  quantities, 

such  as  Nusselt,  Reynolds,  and  Prandtl  numbers.  Such  correlations  for 

single  phase  flow  are  not  adequate  because  these  numbers  are  different 

for  eich  ^hase  and  do  not  take  into  account  the  configuration  of  the 

liquid -vapor  flow.   The  dimensionless  X   is  used  in  an  attempt  to 

tt 

correlate  the  known  properties  of  each  component  with  pressure  drop 
and  heat  transfer  characteristics. 


CHAPTER  II 
TWO  PHASE  FLOW 

A  number  of  studies  have  been  made  of  the  comnlex  configurations 
of  two  phase  flow.  Thece  studies  inviriably  are  broken  down  to  take 
into  account  the  differences  in  the  flow  configuration.  A  designation 
is  given  to  each  type  of  flow.  The  flow  configuration  actually  varies 
in  many  ways  but  several  combinations  of  the  simple  designations  will 
suffice  to  discribe  the  flow  pattern  in  this  investigation.  These  basic 
configurations  are: 

(a)  clear  — —  dry  vapor, 

(b)  mist tiny  water  droplets  carried  along  in  the  main 

vapor  flow  at  low  moisture  contents, 

(c)  annular  flow  —  characterized  by  the  fact  that  the  liquid 

flows  in  an  annulu's  along  the  tube  wall,  \stiile  the 
vapor  passes  at  a  much  higher  velocity  through  the 
center  of  the  tube, 

(d)  slug  flow  —  alternate  slugs  of  liquid  and  vapor  pass 

through  the  tube, 

(e)  bubble  flow  —  small  vapor  bubbles  pass  individually 

through  the  tube  at  low  vapor  contents, 

(f )  pure  liquid  flow  —  no  bubbles. 

These  types  of  flow  were  observed  from  top  to  bottom  as  the  percent 
moisture  wis  increased  from  zero  to  100;".  The  transition  between  two 
types  of  Qow  was  an  over  lapping  affair  and  depended  on   the  flow 
rate  as  well  as  the  moisture  content  (slug  flow  occured  at  lower  percentage 
of  moisture  at  the  lower  flow  rates).  For  example  at  high  moisture 


contents  with  annular  flow  there  would  invariably  pass  an  occasional 
slug  with  rather  consistent  regularity,  and  at  low  moisture  contents 
when  annular  flow  was  starting  to  occur  there  was  certainly  still  a 
definite  mist  in  the  center  of  the  tube.  The  description  of  the  flow 
configuration  definitely  depends  on  the  observer.  The  writer  chose 
to  keep  the  description  simple  and  only  of  secondary  nature  in  this 
investigation.  The  flow  designations  above  were  first  described  by 
Bergelin  (1£)  for  flow  through  a  vertical  tube.  A  considerably 
different  configuration  was  described  by  Martinelli  (9)  for  flow  in  a 
horizontal  tube.  The  basic  difference  for  horizontal  flow  is  that  the 
liquid  stays  on  the  bottom  of  the  tube  while  the  vapor  travels  at  high 
velocity  through  the  top  portion.  Even  so,  at  certain  velocities  an 
assymmetrical  type  of  annular  flow  was  noted  in  the  horizontal  tube. 
This  configuration  difference  accounts  for  the  differences  in  results 
obtained  by  Martinelli  (8)  and  Dengler  (2). 

A  survey  of  the  literature  on  two-phase  flow  has  recently  been 
made  by  Isbin  (U).  The  two  most  widely  used  methods  of  predicting 
two-phase  pressure  drop  are  the  Martinelli  correlations  and  the  friction 
factor  methods.  Of  the  two  methods,  the  Martinelli  correlations  have 
received  the  most  support  and  have  been  chosen  for  comparison  in  this 
investigation.  However,  Isbin  notes  several  drawbacks  to  the  Martinelli 
method.  The  total  flow  rate  parameter  is  not  adequately  provided  for 
and  the  system  pressure  parameter  required  further  treatment. 

A  good  general  discussion  of  the  work  of  some  of  the  recent  develop- 
ments is  included  in  the  latest  edition  of  He  Adams  (8). 


CHAPTER  III 

EQUIPMENT 

The  exnerim  etital  layout  ^resented  in  the  schematic  iagram  of 
Figure  1,  was  obtained  by  modification  of  the  set  up  used  by  Nelson 
(ID. 

The  equipment  consisted  essentially  of  a  flow  system  in  which  steam 
was  taken  from  the  irnin  steam  line,  through  a  gate  valve,  a  small 
centrifix  seperator,  and  a  throttling  valve,  passed  vertically  do-mward 
through  V~>   feet  of  1-|  inch  oiping  and  then  vertically  upward  for  8  feet 
(lfr  inch  oioe),  through  three  feet  of  flexible  hose  (|  inch),  a  sight 
section  (fr  inch),  another  three  feeH  of  flexible  hose  to  the  test  section. 
The  oumose  of  the  flexible  hose  was  to  allow  for  expansion  and  prevent 
undue  stress  on  the  sight  section.  Saturated  water  or  wet  steam  was 
injected  into  the  stream  on  the  lovrer  end  of  the  upward  flow  section. 
Early  in  the  test  a  nozzle  w->s  used  to  spray  the  water  into  the  mnin 
ste^m  line.  Later  it  wis  found  that  a  greater  range  of  moisture  content 
due  to  the  increased  stability  of  the  system  could  be  obtained  v/ithout 
the  nozzle,  with  no  a^oarent  change  in  the  flow  istribution  of  the 
system.  The  tmr^ose  of  the  loop  in  the  steam  line  was  to  allow  the 
steam-wter  mixture  to  approach  equilibrium  in  a  long  straight  section 
before  entering  the  test  section. 

At  the  exit  of  the  test  section  the  flow  was  reversed  and  passed 
vertically  downward  through  a  Centrifix  Type  HA  Separator.  From  the 
steam  outlet  of  trie  centrifix  the  dry  steam  was  passed  thror-h  a  Fj  gher 
and  Porter  Co.,  Series  ^0  Flo-.nr->tor  meter,  then  through  a  needle  control 
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valve.   (This  valve  was  changed  from  preceding  the  steam  flowmeter  ps 
was  used  in  Nelson's  investigation  (11),  to  after  the  flowmeter  in 
or ter  to  insure  dry  saturated  vapor  in  the  flowmeter. ).  After  the  steam 
exit  control  valve,  the  steam  passed  through  16  feet  of  3/U  inch  pipe, 
35>  feet  of  1-  inch  pipe,  seven  feet  of  two  inch  pipe  and  discharged  into 
a  condenser  operated  at  a  slight  vacuum.  From  the  water  outlet  of  the 
centrifix  the  saturated  water  passed  through  a  cooling  coil,  a  needle 
valve,  a  Fisher  and  Porter  Series  700  Flowrator  meter,  and  discharged 
directly  to  the  condenser.  The  needle  valve  ws  adjusted  until  the  water 
level  of  the  centrifix  remained  constant  before  noting  a  reading  on  the 
flowmeter. 

The  moist'jre  content  of  the  steam  leaving  the  Centrifix  was  determined 
by  a  U-Path  Steam  Calorimeter.  The  moisture  content  remained  below 
\%   exceot  at  hi^h  inlet  moisture  contents  with  high  flow,  where  the 
steam  leaving  the  seperator  approached  one  percent  moisture. 

Steam  supply  came  from  a  Babcock  and  Wilcox  F  M  Boiler  at  an 
operating  nreasure  of  200  psig. 

Condenser  condensate  was  pumoed  to  a  suoply  tank  for  a  source  of 
cold  oure  water.  The  water  was  partly  bypassed  through  a  Duolite 
deionizer  on  the  way  to  a  centrifugal  vane  type  pump,  then  through  a 
Fisher  and  Porter  Series  700  Flowrator  meter,  through  a  needle  control 
valve  to  a  heating  section.  The  heating  section  consisted  of  two 
Schutte  and  Koerting  Co.  jet  type  heaters  mounted  in  series.   A  separate 
line  was  run  from  the  main  steam  header  to  sunnly  heating  steam.  The 
heaters  were  more  than  adequate  for  all  conditions  of  these  tests. 
From  the  heater  section  the  saturated  water  (or  under  most  conditions 
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wet  steam)  passed  on  to  the  lower  en  :  of  the  up  r   'low  section  where 
it  was  ::ixed  \dth  the  main  steam  flow  as  mentioned  earlier  in  this  chapter; 

The  tent  section  used  by  Nelson  (11)  was  reworked.   All  new  ther- 
mocouples -md  fittings  were  installed.  Ho.  30  ^age  co^aer-constantan 
thermocouples  were  installed  as  indicated  in  Figure  2.  A  Conax 
Thermocouple  gland  seal  was  used  to  seal  the  stream  thermocouples. 
Pressure  taDS  were  installed  at  the  entrance  and  exit  of  the  test  section 
and  it  three  nlaces  in  the  test  section  as  indicated  in  Figure  2. 

Heat  supply  to  the  sec 'ion  consisted  of  four  independent  heating 
elements  containing  about  70  feet  of  No.  17  Wichrome  V  wire,  wound 
around,  each  of  the  four  sections  in  a  single  layer.  Each  heating  element 
eoimletely  covered  the  section.   A  thin  layer  of  mica  helped  insulate 
the  section  electrically  from  the  heating  wire.  The  power  supply  to 
the  heaters  was  controlled  by  four  two  -gang  variac  assemblies  and 
measured  by  portable  wattmeters. 

Heat  insulation  of  the  terst  section  was  accomplished  by  wrapping 
several  layers  of  glass  tape  on  the  outside  of  the  heating  coils,  cover- 
ing with  a  thick  layer  of  asbestos  cord,  and  over  this  was  fitted  a  2\ 
inch  layer  of  magnesia  brick,  and  a  -|-  inch  of  wet  magnesia  mix.  Cheese- 
cloth was  wrapped  over  the  insulation  and  held  with  cornstarch  paste  to 
prevent  flaking  off.  The  centrifix  and  all  piping  up  to  the  steam  flow- 
meter was  also  insulated.  Heat  loss  was  determined  to  be  negligible 
through  the  insulation. 

A  differential  manometer  was  used  to  measure  pressure  drop  in  the 
section.  The  inlet  pressure  tap  was  connected  to  one  side  of  the  mano- 
meter as  a  reference.  A  pressure  gage  was  also  connected  to  this  tap 
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to  measure  inlet  pressure.  To  the  other  side  of  the  manometer  was  con- 
nected a  manifold  to  which  the  other  four  pressure  taps  were  connected. 
The  manometer  lines  were  connected  at  the  test  section  by  llorgren  needle 
valves  and  at  the  manifold  by  quick  acting  toggle  valves  manufactured  by 
Hoke.  During  the  experiment  the  pressure  tap  lines  were  kept  full  of 
water  to  insure  a  known  water  level.  To  keep  these  lines  full  of  water, 
wet  rags  were  applied  to  the  lines  to  condense  any  steam  entering. 

It  has  been  determined  as  a  result  of  this  experiment  that  it  would 
have  been  better  to  use  the  center  pressure  tap  as  the  reference  pressure. 
Because  of  end  leakage  of  heat,  only  the  second  and  third  sections  were 
used  for  heat  transfer  calculations.  At  high  flow  rates  the  pressure 
drop  from  the  inlet  to  the  center  is  significant.  By  moving  the  inlet 
pressure  gage  to  the  center  of  the  test  section,  the  average  test  section 
pressure  can  more  easily  be  held  constant. 

All  components  of  the  system  were  constructed,  as  far  as  practicable, 
of  non-ferrous  material. 
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CHAPTER  IV 
OPERATING  PROC  ID  E 

At  the  beginlng  of  each  operating  oeriod  the  test  section  was 
filled  with  cold  vnter  from  the  supoly  tank  and  the  manometer  system 
ws  filled  with  water  by  bleeding  off  any  bubbles  in  the  line,  accumulated 
since  the  previous  operating  oeriod.  Then  the  system  was  flushed  by 
flowing  a  large  quantity  of  steam  through  the  test  section  to  the  con- 
denser. The  con "ens ate  was  allowed  to  run  down  the  drain  until  it  was 
clear  an'  free  of  particles.  The  drain  was  then  secured  and  the  con- 
densate by-passed  to  the  storage  t~>nk.  This  procedure  required  about 
two  hours. 

The  operating  variables  that  could  be  controlled  were  pressure, 
heat  flux,  E*low  rate,  and  moisture  content.   As  noted  in  Appendix  II 
certain  limitations  of  flow  existed  Tor  each  oressure  held  at  the  test 
section.  The  inlet  pressure  to  the  test  section  was  maintained  at  10U 
psig  for  all  runs  except  three  runs  at  Hi  9  psig.   The  heat  inout  oer 
section  was  limited  to  1.5>  Kl>T  for  all  runs.   At  one  KV7  per  section  t   e 
experimental  error  was  magnified.   At  two  KW  per  section  the  circuit 
drew  excessive  current, 

A  totil  of  69  runs  were  made  of  which  ^0  were  considered  to  be 
satisfactory,  k   summary  of  the  da1   Ls  listed  in  Appendix  I.  Unsatis- 
factory.runs  were  cai  sed  by  instability  of  the  system  making  it  difficult 
to  maintain  equilibrium  con1! "ions  at  a  desired  oressure.  The  unsati  - 
factory  runs  are  not  inclu  led  Ln  th<     n  ix. 

The  desired  ["low  rate  and  moisture  content  was  principally  obi 
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equilibrium.   For  the   at  urn<      -   KW 

oer  section  n  'ter  the  liqui'1  'low  h    een  est'  lj    .  eric 

of  time,  rluri  .'  straent:-  were  n-  ,  the  h rere 

increase  '  to  one  KV7  ner  s^c'ion  an n  then  1-|  KV7  per  section.  If  the 
increase  was  too  raoid,  excessive  current  caused  the  fuses  to  blow. 
The  time  required  to  adjust  to  5  esire  ooint  varie ■'   'om  30  minutes 

i"or  on  isothermal  run  up  to  on  hour  or  ;.iore  .  r  ah    d  run.  The 
stability  of  the  ■  stem  ^r<  ttly  affected  the  time  re   .   to  adjust  to 
a  'esire'"  value.   n:e  system  approached  instabili1     low  'low  rites 
and  at  hi  ;h  moisture  content  runs  for  all  flows  as  the  system  approac 
'  lug  "l^w.  The  stability  was     tly  Improve   over  the  whole 
when  the  ---tor  nozzle  was  removed.  The  stabilit;   or  a  desired  co  n  ition 
could  soro  times  be  imoroved  by  reducing  the  heal  Lng  liquid  to  below 
saturation  and  replacing  it  with  ste?yn  passing  through  t'     i   steam 
line.   It  other  bim<  s  the  main  steam  lim   -•-  •  out  off  al together j  the 
er  was  he^ti  d   Lth  one  heat<  r,  ti   the  oressur*  w  s  ■   ntrolled  by 
the  throttle  valve  to  the  other  heater. 

For  each  run,  thermocouple  millivolt  readings,  test  section  i  let 
pressure,  ifferential  pressure  droos,  inter  lowmeter  re     ,         am 

'lowm     *ea  ling,  and  pressure,  and  heat  inout  in  KW  were  reco 

Several  additional  n  s  were  m     i     to  check  "■  ,  lidit 
data.  Tie  steam  Flow  met  c       s  c  librate  .  The  • 

tained  in  t  o  Inches  of  the  centrifix  w-.s  1         or  use  at  very 
low  moisture  content,  nd  exten:  Lv  :omparisons  of  thermc  Lngs 

■   nde.   Pure  liq"i  ■  <     to  c\    luai    n  1     t 

validity  of  heat  transfer  calculations.  These  calibration  runs  are 
more  thoroughly  described  in  Appendix  II. 
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CHAPTER  V 
METHOD  0^  CALCULATION 

The  film  heat  thransfer  coeficient,  h,  is  defined  as  the  pro- 
portionality factor  in  Nekton's  Law  of  Cooling  (8): 

dq  :  h  dA  (tg  -  t) 
The  surface  temperature  (t  )  was  determined  from  the  average  of  the 
thermocouple  readings  in  the  second  md  third  sections,  talcing  into 
account  the  temperature  difference  between  the  surface  pjid  the  location 
of  the  thermocouples.  This  temperature  difference  was  determined  to 
be  5.02  °F,  for  Q/A  equal  to  7% 300  BTU/hr-ft2  (corresponds  to  1.5 
KVJ  oer  section),  from  the  Fourier  conduction  equation  for  a  circulnr 

tube  (8): 

2xk  L  (ti  -  to) 

q  z 


In  r^ri 

The  temoerature  of  the  fluid  (t)  in  Norton's  equation  was  taken  as 
the  saturation  temperature  corresponding  to  the  pressure  at  the  center 
of  the  test  section.  Because  of  thermal  entrance  effects  on  the  first 
section  and  end  leakage  from  the  first  and  last  sections,  only  the  two 
center  sections  were  used  in  the  final  calculations. 

The  pressure  at  the  center  of  the  test  section  was  determined  from 
the  inlet  pressure  gage  reading  les;  the  pressure  drop  indicated  by 
the  differential  manometer. 

The  moisture  content  was  taken  at  the  exit  of  the  test  section. 
The  amount  of  liquid  evaporated  (about  20  lbs/hr)  was  very  small  in 
proportion  to  the  total  flow. 

The  single  phase  points,  pure  liquid  and  pure  v  por,  were  calcu- 
li ted  from  the  convection  equation  for  "low  in  tubes  (8) I 
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(  hD/k  )  =  0.023  (nc),8(Pr)Jl 
The  calculations  for  the  liquid  phase  (hr)  were  use   or  the  comparison 
with  Dengler's  results  (h/h^  vs  l/X^)*  Fibres  k   and  5. 

The  pressure  drop  was  measured  by  a  differential  manometer  across 
the  second  and  third  sections.  The  single  phase  presrure  drops  were 
calculated  from  standard  friction  factor  charts. 

The  dimensionless  parameter  X^  was  developed  from  dimensional 
analysis  by  Ilartinelli  (9): 

W/s0.9  ,  f*,x0.S  ,^ox0.1 


h*  •  (^)°-9  (^)°-5  ^ 


It  has  been  found  useful  in  correlating  data  of  tv;o  phase  flow  over 
a  wide  variety  of  conditions.  X^  was  calculate."  as  a  function  of 
X  for  two  different  pressures  (lOh  psig  and  Ui9  psig)  with  data  from 
the  steam  tables  for  correlation  of  the  data  of  this  investigation. 

It  appeared  to  be  possible  to  predict  values  of  h/h..  from  measured 
values  of  R^,  or  v^sa   versa,  with  the  assumptions  that  there  is  annular 
flow,  that  there  is  no  nucleate  boiling,  and  that  an  equivalent  pipe 
size  can  reqresent  the  cross  section  of  the  liquid  area.  An  equivalent 
diameter  De  is  defined  to  represent  the  equivalent  pipe.  For  single 
phase  forced  convection:   (using  equivalent  pipe). 

hDeA  =  0.023  (Ree)'8  (Pr),U 

where  Ree  =  De  m  (l-x)/A/*-, 

or  to  obtain  h/hL  rearrange  and  divide  by  h  in  the  same  form: 

h/hL  =   0.023  (k/Dfi)  (ReP),J  (Pr)* 
0.023  (k/D)   (Re)'3  (Pr)'*4- 
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Simplifying, 


D2-2g(l-x) 


Qe 


o  or 


.8 


Since  DQ  =  D  |R^  by  definition, 

.3 


h/hL  = 
or  rearranging, 


(1-*) 

TTfe 


.889 


Similar  attempts  to  calculate  pressure  drop  were  not  successful. 
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CHAPTER  VI 

RESULTS  AND  CONCLUSIONS 

The  measured  he't  transfer  coefficient  vs  oercont  moisture  by 
weight  is  shown  in  Figure  3.   It  is  noted  that  in  the  annular  flow 
region  (about  30$  to  80'b  moisture)  the  variation  is  essentially 
linear*   From  the  initial  observation  of  these  results  it  appears  that 
a  rclitively  unpredictable  mechanism  is  causing  an  increase  in  the  local 
film  heat  transfer  coefficient.   But  the  results  of  these  tests  combined 
with  the  vrork  of  Dengler  (2)  indicate  that  the  local  film  heat  transfer 
coefficient  for' vertical  two  phase  flow  is  still  a  single  phase  phenomena 
in  the  annular  flow  region.  The  only  significant  resistance  to  heat 

"low  is  the  liquid  laminar  boundary  layer.  For  the  heat  flux  used  in 
this  investigation  no  nucleate  type  of  boiling  occured  at  the  metal 
surface,  and  the  substantial  increase  in  heat  transfer  coefficient 
can  be  attributed  entirely  to  the  increase  in  liquid  velocity  and  the 
consequent  reduction  of  the  laminar  boundary  layer  thickness. 

It  will  also  be  noted  in  Figure  3  that  the  presage  is  a  very 
significant  variable.  Dingier  (2)  has  suggested  plotting  h/h^  vs 
1/X^  and  this  brings  together  the  runs  at  different  flow  rates  (Figure 
h)   as  well  as  at  different  pressures.  The  pressure,  however,  h  s  riot 
been  completely  provided  for  as  is  seen  by  comparing  the  results  (Figure 
J4)  of  this  investigation  at  lOU  psig  with  Denver's  results  at  -bout 
atmospheric  pressure.  Difference  in  oioe  size  (Dengler  used  1"  diameter 
comnared  to  i"  diameter  for  this  invesl    ,.on)  may  partly  cause  the 

inference  in  results. 
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Dengler  (2)  also  measure  (R^)  the  volumetric  friction  of  liquid 
(See  Chapter  I).  Using  Dengler's  values  of  R-r  (Figure  £)  and  a  pressure 
of  Hi. 7  psia,  h/hr  was  calculated  from  single  nhnse  considerations 
(See  Chanter  IV  for  method  of  calculation)  and  plotted  on  Figure  U, 
The  deviation  can  be  attributed  mainly  to  the  degree  to  iihich  an  equiva- 
lent pipe  oes  not  represent  the  actual  flow  at  the  surface  of  the  tube. 
In  Dealer's  experiments  the  oressure  varied  from  7  psia  to  29  psia  and 
this  may  al  o  cause  some  of  the  discrepancy.  Even  with  the  discrepancy 
as  shown  it  is  felt  that  the  hypothesis  that  the  film  he  it  transfer 
coefficient  is  n   single  phase  problem  is  substantiated. 

To  carry  the  comparison  one   step  further,  Rt  was  calculated  from 
the  experimental  points  (h/hy)  of  this  investigation  and  olotte  on 
Figure  5,  A  definite  shift  upward  is  noted  for  the  increase  in  oressure 
This  trend  is  also  noted  by  iiartinelli  (11)  >r  horizontal  flow,  but 
only  in  a  qualitative  manner.   This  inv  stigation  suggests  th-^t  the 
shift  in  Rt  (also  h/hT )  is  much  greater  in  the  lox.  oressure  regions 
than  at  higher  pressures.  Ibis  further  suggests  that  a  modification 
of  the  density  term  in  X. ,  is  required  to  give  a  better  corr-     on. 

The  measure1  values  of  pressure  droo  nor  unit  length  vs  percent 
moisture  are  shown  in  Figure  6.  This  pressure  droo  is  composed  o 
three  terms : 

1.  The  frictlonal  two  phase  pressure  drop,  (AP/L)Tpp. 

2.  The  pressure  drop  due  to  change  in  elevation,  "az"/L, (psi/ft). 

3.  The  pressure  'roo  hie  to  change  of  momentum  during  evaporation. 
This  value  was  negligible  in  V  is  investiq-'tion  because  o!   the  low 
vaporization.  This  fact  is  varified  by  the  indistinguishable 
difference  between  the  heated  ^r\ti   the  isothermal  ^oints  '    :*e  6), 
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For  correlation  of  pressure  drop  Martinelli  (7)  ^n  (11)  has 

sted  the  use  of  .1  trimeter  0^       where:        <$2         _     (AP/l)"ppp       . 

Ut  ~  (4P/L)t 

The  numerator  represents  the  friction-il  component  of  two  phase  flow  and 

the  denominator  represents  the  friction^!  pressure  ;roo  Cor  an  equal 
mass  rate  oP  liquid  ilonc.   In  Figure  7t  fi^       vs  X^  is  plotted  for 
this  investigation  and  compared  with  the  results  of  Den^ler.  The  dashed 
line  indicates  the  nossible  range  of  Dengler's  c\irve  at  the  left  en  . 
The  difference  of  results  is  more  oronounced  in  this  comparison  than  is 
easily  explained  by  the  difference  in  onerating  pressure.  Mar  L  elli 
(11)  ag''in  shows  a  trend  in  this  direction  with  increasing  oressurc  but 

not  to  this  extent.   It  is  believed  by  the  writer  that  the  pioe  rize 

i 
becomes  more  important  for  the  pressure  drop  correlations.  This  added 

to  the  difference  in  onerating  ores sure  may  account  for  the  large 

discrepancy. 

Initial  attempts  to  calculate  friction al  ores sure  drop  from  single 

ohnse  considerations  (knowing  Rt)  were  unsuccessful.  It  is  believed 

that  agnin  oi^e  size  is  the  significant  variable  that  was  not  adequately 

t^ken  into  account. 

Summary: 

1.  The  film  heat  transfer  coefficient  is  linear  in  the  annular 
flow  region  when  plotted  against  percent  moisture. 

2.  The  only  significant  resistance  to  neat  flow  is  the  "iaiid 
laminar  boundary  layer,  and  hence  the  increase  in  h  is  due  primarily 
to  liquid  velocity  increase. 

3.  The  ratio  h/hj  vs  lAtt  correlates  well  and  comoares  favor 
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vrlth  results  of  Dengler  (2). 

h.     The  ntio  h/hy  can  be  calculated  approximately  from  PL  data 

(in  the  annular  flov;  region)  by  modifying  single  phase  methods. 

5.  Effects  of  oressure  on  h/hT  are  much  greater  for  the  low 

L       ~4 
pressure  regions. 

2 

6.  The  pressure  dror>  narameter  0-i   correlates  well  with  X.  .  . 
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CHAPTER  VII 
RE(     (DA^TONS 

Several  significant  problems  are  suggested  by  the  results  of  this 
investigation: 

1.  The  Prediction  of  two-phase  film  heat  tranfer  coefficients 
from  the  direct  measurement  of  the  volume  fraction  of  liquid  in 
an  isothermal  experiment, 

2.  Conversely,  the  indirect  measurement  of  the  volume  fraction 
of  liquid  (particularly  at  higher  pressures)  from  the  measurement 
of  the  loci  film  heat  transfer  coefficient.  Along  with  this 
investigation  the  effects  of  smaller  test  section  'iameter  could 
be  investigated. 

3.  The  analytical  prediction  of  two  Phase  ores sure  droo  from  a 
single  phase  approach.  This  would  assume  that  the  frictional 
pressure  drop  is  due  to  the  increase  in  the  liquid  velocity  onl;;. 
Certain  changes  are  suggested  in  the  construction  of  the  test 

section  to  avoid  some  of  the  difficulties  encountered  in  this  investi- 
gation. The  reference  pressure  and  the  inlet  pressure  gage  should 
be  moved  to  the  center  of  the  test  section  since  this  is  the  significant 
pressure  in  the  ultimate  results.   The  inlet  and  outlet  pressure  taps 
should  be  moved,  from  the  vicinity  of  the  stream  thermocouples  to  obtain 
consistent,  readings  of  Pressure  differences  alon  -  the  test  section. 
It  is  recommended  that  these  pressure  taps  be  installed  immediately 
before  and  after  the  heated  sections  giving  equal  spacing  between  all 
pressure  taPSt  Difficulty  was  encountered  with  leakage  at  the  pressure 


27 


tro   connections.  This  was  caused  at  the  point  inhere  the  hard  br?.ss 
fitting  into  the  soft  test  section.  T'e  thread  could  be  strip  c   by 
installing  the  fitting  only  "inger  tight.  It  is  recommended  that  more 
threads  and  a  slightly  larger  diameter  of  fitting- be  used.  It  shovld 
be  remembered  vrhen  lesigning  the  test  section  that  the  first  and  last 
sections  -rill  not  p;ive  significant  results  because  of  end  leakage  and 
a  thermal  entrance  ef:'ect. 
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APPENDIX  II 
EQUIPMENT  LIMITATIONS 

Flow  'l?te  limitations 

It  w^s  found  th-'t  certain  equipment  limitations  prevented  the 
OPern ' ion  at  certain  desired  flow  rates  and  moisture  contents.  Early 
in  the  investigation  it  was  necessary  to  ascertain  the  limitations  of 
:.■  equipment  in  order  thai  the  widest  Possible  range  of  data  coull  be 
taken.  Figure  8  shows  tl  e  upper  limit  of  flow  for  various  test  section 
inlet  pressures.  These  were  obtained  by  varying  the  amount  of  saturated 
liqui  1  introduced  with  the  main  steam  inlet  valve  wide  open.  One 
additional  curve  could  have  been  obtaine  I  at  100  psig  but  the  steam 
exit  throttle  valve  lost  control  in  this  range,  with  zero  moisture. 
For  certain  moisture  contents,  data  could  be  taken  down  to  on  inlet 
pressure  of  75  ^sig.  The  lower  portion  of  the  curve  which  is  cross 
hatched  represents  the  lower  limit  on  the  scale  of  the  steam  flowmeter. 

The  approximate  points  obtained  in  this  investigation  are  indicated. 
For  700  lbs/hr.,  data  were  also  taken  for  three  points  at  UU9  psig 
inlet  pressure. 

Steam  Flowmeter 

The  "Factor  Tag"  indicated  100$  reading  on  the  steam  flowmeter  was 
equivalent  to  11^0  lbs/hr  at  20(3  psig.   For  other  values  of  pressure 
the  full  scale  reading  had  to  be  modified  as  a  function  of  density. 
The  following  formula  was  used  to  determine  the  actual  flow: 


W  =  (Scale  Reading)  (11^0)  \/_£_ 
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The  d  nsity  was  obtained  from  the  specific  volume  data  in  the  steam 
tables.  The  metering  tolerance  was  estimated  by  the  manufacturer  to 
be  about  two  to  three  percent.  The  scle  could  be  read  within  a  quarter 
of  a  division  ('"). 

A  calibration  check  was  made  at  l'i9  Psig  by  holding  the  steam  flow 
constant  and  measuring  the  condensate  for  a  period  of  time  in  a.  weight 
tank. 

'Jatcr  Flowmeter 

The  water  flowmeter  float  was  designed  for  2.7  GPM  at  full  scale 
(100%').  This  is  equivalent  to  1350  lbs/hr.  The  liquid  was  cooled  to 
about  60°F  before  passing  through  the  flowmeter  and  therefore  no  temp- 
erature correction  was  required.  The  liquid  "low  meter  could  not  be 
used  Cor  low  flow  rates.  To  measure  low  flow  rates  the  capacity  of  two 
inches  of  the  centrifix  sight  glass  was  determined  to  be  1.55  rounds. 
For  low  flow  rates  time  required  for  the  liquid  to  rise  two  inches  in 
the  centrifix  sight  glass  (with  the  liquid  exit  throttle  valve  secured) 
was  measured  and  the  flow  rate  thereby  determined. 

Manometer  Correction  Factor 

The  manometer  readings  were  measured  in  inches  of  fluid  with  an 
equivalent  specific  gravity  of  12.6  (specific  gravity  of  mercury  minus 
the  soecific  gravity  of  water )#  Because  of  the  difference  in  the  liquid 
level  of  the  different  oressure  taps  a  further  correction  factor  had  to 
be  added.  This  was  taken  as  the  measured  distance  (in  inches)  from  the 
inlet  pressure  tap  to  oressure  tap  concerned,  divided  by  12,6.  The 
calculated  correction  factor  was  verified  by  taking  readings  with  zero 
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flow.  In  interpreting  the  data  all  pressure  readings  were  plotted  to 
determine  consistency.  If  the  readings  were  consistent  the  corrected 
pressure  drop  across  the  tvro  center  sections  was  divided  by  1.033  feet 
and  multiplied  by  .U5>5  to  obtain  (A.P/L)  in  p si/ft.  It  was  found  early 
in  the  test  that  there  was  inconsistency  in  the  readings  at  the  inlet 
pressure  tap  and  the  outlet  pressure  tap  due  to  the  presence  of  the 
stream  thermocouples  at  the  same  location.  Later  in  the  tests  the  stream 
thermocounles  were  removed  and  consistent  results  were  obtained. 

Thermocouple  Readings 

Numerous  thermocouple  readings  were  ^lotted  for  both  isothermal 
and  the  heated  runs  and  compared  with  the  saturation  temperature 
corresponding  to  the  pressure  at  each  point  along  the  test  section. 
While  there  was  some  scatter  in  these  results  it  was  felt  that  the 
average  of  the  thermocouple  readings  of  the  two  center  sections  gave 
a  reliable  value. 
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APPENDIX  III 
ADDITIONAL  CALCULATIONS 

Single  Phase  Flow 

Attempts  to  make  heated  single  phase  runs  were  unsuccessful.  For 
the  liquid  run  there  was  insufficient  nixing  to  measure  the  exit  teraoer- 
ature  of  the  liquid.  Thermal  entrance  effects  were  greatly  magnified. 
End  leakage  was  evaluated  and  it  was  confirmed  that  the  data  in  the  two 
end  sections  was  unreliable.  For  the  vapor  run  the  fluid  was  super- 
he  ote  >  in  the  first  section  which  did  not  give  sufficient  data  to 
calculate  reliable  values  of  h.  Therefore  the  single  phase  values 
were  calculated  and  are  listed  below. 


Liquid 

Mass  Inlet 

Rate  Pressure 

(lbs/hr)  (r>sig) 


liOO 

700 

1000 

1250 

700 


loii 


111? 


hL   (AP/L),    ("iZ»/L)i 
(psi/ft)  (psi/ft) 


872  .0078 

1338  .0216 

1812  .Oh 12 

2170  .0628 


•U33 
.U33 


1381;  .0215    .0215 


Vapor 


280 
U38 
583 
969 

U33 


(*P/L)g 

(psi/ft) 

1.165 
3.37 
6.1|8 
9-U8 

2.U7 


(«4z»/L) 


g 


0 
0 
0 
0 


Thermal  conductivity  of  copper 

The  coefficient  of  thermal  conductivity  was  obtained  from  the 
manufacturer,  k  equal  to  222  BTU/hr-ft-°F.  This  value  was  extrapolated 
from  a  comparison  of  the  variation  of  k  for  pure  copper  with  temperature 
to  k  equals  215  BTU/hr-ft-°F  at  350°F.  This  change  had  little  effect 
on  the  results. 


/ 
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Insulation  Lo3S 

The  insulation  loss  was  calculated  as  follows! 

Let  AT  =  70  °F  (Temperature  of  insulation  about  1£0  °F) 
For  Vertical  plates  (L>1  ft)i  h  =  .3(4T)*2  '  "  »3(70)#2 

h  =  0.87  BTU/hr.ft2-°F 
The  heat  loss:  q  =  h  A  AT  "  (.87)(*/8)  (70) 

q  =  23.9  BTU/hr/section  ***     .007  KW 
The  temperature  of  the  insulation  did  not  reach  this  value,  therefore 
the  insulation  loss  was  assumed  to  be  negligible. 
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